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Abstract. The coherent torques between stars on orbits near massive black holes (MBHs) lead 
00 ' to resonant angular momentum relaxation. Due to the fact that orbits are Keplerian to good 

approximation, the torques efficiently change the magnitude of the angular momenta and rotate 
the orbital inclinations. As a result the stars are rapidly randomized. The galactic MBH is a 
good system for the observational study of resonant relaxation. The age of the young B-stars at 
a distance of ~ 0.01 pc from the MBH is comparable to the resonant relaxation time, implying 
that resonant relaxation may have played an important role in their dynamical structure. In 
contrast, the O-stars in the stellar disks at ~ 0.1 pc are younger than the resonant relaxation 
time, as required by their dynamical coherence. Resonant relaxation dynamics dominates the 
event rate of gravitational wave (GW) emission from inspiraling stars into MBHs of masses 
' comparable to the Galactic MBH. Resonant relaxation leads to rates < 10 times higher than 

f"^ , those predicted by 2-body relaxation, which would improve the prospects of detecting these 

events by future GW detectors, such as LISA. 
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52 ' 1. Introduction 

In "collisional" stellar dynamics the potential in which a star moves is considered to be smooth 
to first order, and the fact that the potential is in fact made out of discrete stars is treated as a 
^ ■ perturbation (e.g. Chandrasekhar Binney &i Tremaine ^0])- Stars have orbital parameters 

such as angular momentum and energy that are conserved in the smooth potential, and these 
parameters indeed remain constant for many dynamical times td- Only the weak encounters 
allow for changes in these quantities. For example, two stars which interact with each other, 
can exchange energy and angular momentum, so that after their encounters their orbits are 
described by slightly different quantities. In this way not only the orbits of individual stars 
are modified, but the distribution function (DF) of the entire system can change, and evolve 
towards a steady state. The time-scale over which a system evolves is the relaxation time t r . In 
most systems td <C t r , and relaxation can indeed be treated as a second order effect. 

Analyses of the evolution of the DF near a MBH have almost exclusively relied on the 
assumption that the mechanism through which stars exchange angular momentum and energy 
is dominated by uncorrelated two-body interactions 1 . Any encounter is assumed to be unrelated 
to previous and future encounters, and changes in energy and angular momenta are considered 
to be drawn from a specified random distribution. Relaxation can therefore in a meaningful way 

1 iV-body simulations form an exception; for iV-body simulations near MBHs see e.g. Baumgardt, Makino & 
Ebisuzaki 0,0; Preto, Merritt & Spurzem g3|; Merritt & Szell (S3 



be considered to be a random walk process. In the context of stellar dynamics near MBHs, the 
assumption of uncorrelated encounters is made in Fokker-Planck models (e.g. Bahcall & Wolf 0; 
Bahcall & Wolf [£] ; Cohn & Kulsrud ; Murphy, Cohn & Durisen [HH] ) , where the microscopic 
interactions are expressed by the diffusion coefficients, and in Monte Carlo simulations (e.g. 
Shapiro & Marchant [12]; Marchant & Shapiro |35j . |36j : Freitag &; Benz |15| . |16j). Stars 
around MBHs are described as moving in the smooth average potential of the MBH and the 
stars, and the scattering by the fluctuating part of the potential is modeled as a hyperbolic 
Keplerian interaction between a passing star and a test star. 

The (non-resonant) relaxation time Tnr can be defined as the time it takes for the energy £ 
of a typical star to change by order unity. This is also the time it takes for its specific angular 
momentum J to change by an amount of order J c (£), the maximal angular momentum for that 
energy. On Keplerian orbits J c = \/GM 9 a, where o is the semi-major axis. The "non-resonant" 
relaxation time Tnr of stars of mass M* can be written in the Keplerian regime as 

r ™ = ^(t) 2 ^) (M - >>M *>' (1) 

where P = 2ir ^ a? / (GM,) is the orbital period and A\ is a dimensionless constant which includes 
the Coulomb logarithm. For some stellar systems (like galaxies), Tnr is much longer than the 
age of the system, implying that the system cannot evolve significantly towards steady state by 
two-body interactions. For other systems, like very dense stellar clusters, the relaxation time can 
be as small as a few Myr, and such systems may even evaporate within a Hubble time. In our 
GC, the relaxation time is somewhat smaller than the age of the system, Tnr ~ fewGyr (e.g. 
Alexander , [I] ) , implying that the system has evolved considerably, and two-body relaxation 
effects such as mass-segregation have occurred (Bahcall & Wolf Freitag, Amaro-Seoane & 
Kalogera ,19, ; Hopman & Alexander 28 ). At the same time, the relaxation time is much longer 
than some other relevant times in the GC, in particular the age of the youngest stars. 

The assumption of uncorrelated two-body interactions is well-justified in many systems, such 
as globular clusters, where stellar orbits are not closed. However, the special symmetry of a 
Keplerian potential leads to closed, elliptical orbits. The fact that the orbits are closed can 
be exploited in numerical treatment, and also leads to unique dynamical features (see also 
the contribution of Touma in this volume). Since t r 3> td, the orbits remain closed for many 
dynamical times, the system may be thought of as a set of "wires" with the mass of the star 
smeared out over the orbits. In this picture, it is the wires that interact and cause the evolution 
of the system, rather than point particles interacting at given locations. The idea is reminiscent 
of the Kozai mechanism in triple stars (Kozai |H()|). 

Rauch & Tremaine |32] first used this approach in the context of many body stellar dynamics 
near MBHs, and coined the term resonant relaxation (RR), after the 1 : 1 resonance between 
the radial and azimuthal frequencies in a Keplerian potential. The wire approximation is only 
relevant for times <C t u , where t u is the time for the orbit to precess. Precession may be caused 
by the fact that the potential is not entirely determined by a point mass, and there is still some 
extended component to the potential due to the stellar mass; this is especially the case far away 
(> 0.1 pc) from the MBH. Closer to the MBH (< 0.01 pc), precession may be dominated by 
effects of General Relativity 2 . 

1.1. Scalar resonant relaxation 

Scalar relaxation results in changes in both the direction and the magnitude of the angular 
momenta. The RR time Trr is estimated by evaluating A.J U , the coherent change in the 

2 For the parameters of interest here, Lense-Thirring precession is much less efficient than mass and GR precession 
even for a maximally spinning MBH, and we do not consider it here. 



magnitude of the specific angular momentum up to a time t u . The change AJ^ is then the step 
size for the non-coherent growth of the angular momentum over times t>t UJ . Two nearby stars 
with semi-major axes a exert a mutual specific torque ~GM 1 /a. To zeroth order the torques 
of the stars on a test wire cancel, so that within a distance a from the MBH the net torque on 
a test star is determined by the Poissonian excess torque J~ \J N(< a)GM+/a and 



AJ W ~ jtu = yjN{<a){GM*/a)t u 



(2) 



For t > t w the torques on a particular star-wire become random, and the change in angular 
momentum grows in a random walk fashion with a timescale Trr~ (J c / AJ U ) t u , defined as 



T RR = A 



RR" 



N(>£)P 2 (£) A RR 



M m 



(3) 



> 2 (>£) t u N(<a) 

where \i = NM±/(M, + NM±), A RR is a numerical factor of order unity, to be determined by 
simulations, and the last approximate equality holds in the Keplerian regime. 

Over most of the relevant phase space the precession is due to the deviations from pure 
Keplerian motion caused by the potential of the extended stellar cluster. This occurs on a 
timescale t u = tu = [M,/N(< a)MJP(a), assuming N(< a)M* <C M,. The J-averaged RR 
timescale can then be written as 



-RR 



A RR ^-P(a) 



A RR M* 



N{<a)T NR . 



(4) 



A A M. 

Since <C Inr for small a where N(< a)M+ <C M,, the RR rate of angular momentum 
relaxation is much higher than the rate of energy relaxation in the resonant regime. This 
qualitative analysis has been verified by detailed numerical iV-body simulations by Rauch & 
Tremaine [IS] and by Rauch & Ingalls |46j . 

For most of parameter space, orbital precession is dominated by the mass of the stellar 
cluster and the RR timescale is well approximated by T RR ~ T RR - However, very close to the 
MBH, or on wide orbits with very low angular momentum, so that the periapse is close to 
the Schwarzschild radius of the MBH, precession is dominated by GR effects. In this case the 
timescale for precession is given by t w = icR = (8/3) ( J/ «/lso) 2 -P; here Jlso = (4GM./c) is 
the angular momentum of the last stable orbit (LSO). When tcR "C £m and GR precession 
dominates, the RR timescale is (Eq. |1J) 
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1 RR 



8 RR \Mj \ J J N(<a)- 



(5) 



Generally, GR precession and mass precession occur simultaneously, and the scalar RR 
timescale T^ R (£, J) is given by substituting l/t u = \l/tM — 1/*gr| i n Eq. ©, where the opposite 
signs reflect the fact that mass precession is retrograde whereas GR precession is prograde. Thus, 
the scalar RR timescale is 
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We use the relation d(J 2 )/J c 2 = dt/T£ R (£, J) (Eqs. © to define the J -averaged time it takes 
a star to random-walk from J = J c (£) to the loss-cone J = J\ c as 



J lc 



dJ 2 T RR (£,J). 



(7) 



1.2. Vector resonant relaxation 

For time scales much larger than the dynamical time, orbits precess and describe a rosette shape. 
One can then consider the torques between different rosettes rather than between different 
wires. Since the rosettes describe planar rings to good approximation, they cannot modify the 
magnitude of the angular momentum of the star, but they can change the direction of the 
angular momentum vector. This process is known as "vector resonant relaxation" (Rauch &; 
Tremaine Vector RR grows coherently {pat) on timescales t<^t v , where t v is the timescale 

for a change of order unity in the total gravitational potential ip caused by the changes in the 
stellar potential (p* due to the realignment of the stars as they rotate by tt on their orbit, 



** = J r-— -- oT-FTTTT^ (8) 



the last approximate equality holds for NM+ <C M 9 . In analogy to scalar RR (Eq. |2J), the 
maximal coherent change in J is |AJ<J ~ Jt v ~ J c , that is, J rotates by an angle 0(1) already 
at the coherent phase. On timescales 4 3>t«, |AJ„| cannot grow larger, as it already reached 
its maximal possible value, but the orbital inclination angle is continuously randomized non- 
coherently (oct 1 / 2 ) on the vector RR timescale (Eq. I3J), 

rpv 0A v NV\>£) pfe ,. n , v (M.\ P(a) 

T RR = 2A RR -77^ P(£) - 2A RR TT \ ' ( 9 ) 



fi(S) \M*J N l l 2 (<a) 

where the last approximate equality holds for NM+ <C M, . 

It is that while the torques driving scalar and vector resonant relaxation are the same, vector 
RR is much more efficient than scalar RR, T^ R <CT^ R , due to the much longer coherence time 
N 1 / 2 tM^ > tM ■ Furthermore, vector RR proceeds irrespective of any precession mechanisms 
that limit the efficiency of scalar resonant relaxation. 

2. The origin of the young stellar population in the Galactic center 
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Figure 1. Stellar components, timescales 
and distance scales in the GC. The NR 
timescale Tnr (top straight line); the 
timescale estimated for 1 M stars (top 

curved line) and 10 M Q stars (bottom curved 
line); the timescale T RR (bottom straight 
line); the position and estimated age of the 
young stellar rings in the GC (filled rectan- 
gle in the bottom right); the position and age 
of the S-stars if they were born with the disks 
(empty rectangle in the bottom left); the posi- 
tion and maximal lifespan of the S-stars (filled 
rectangle in the middle left). Reprinted with 
permission from the Astrophysical Journal 



Figure Q compares the distance scales and the ages or lifespans of the various dynamical 
structures and components in the inner pc of the GC with the relaxation timescales. The NR 
timescale in the GC, which is roughly independent of radius, is TNR~fewx 10 9 yr (Eq. The 
scalar RR time Tj| R is shown for M* = 1, 10 Mq. At large radii the RR time decreases towards 
the center, but for small radii, where GR precession becomes significant, it increases again. The 



vector RR timescale T RR , in contrast, decreases unquenched with decreasing radius. Structures 
with estimated ages exceeding these relaxation timescales must be relaxed. 

Two distinct young stellar populations exists in the GC. At distances of 0.04-0.5 pc from 
the MBH there are about ~ 70 young massive OB stars (M* 3> 10 M , lifespan of t* = 6 ± 2 
Myr), which are distributed in two nearly perpendicular, tangentially rotating disks (Levin &; 
Belobodorov 32 : Genzel et al. 20 : Paumard et al. [E]). It appears that these stars were 
formed by the fragmentation of gas disks (Levin & Belobodorov [32]; Levin jHJ; Nayakshin & 
Cuadra jinji Nayakshin & Sunyaev |41j : Nayakshin |42p. Inside the inner 0.04 pc the population 
changes. There is no evidence for old stars, and the young stars there (the "S-stars") are main- 
sequence B-stars (M*< 15 M Q , lifespans of 10 7 < t*<2x 10 8 yr; Ghez et al. [2T]; Eisenhauer et 
al. J3]) on randomly oriented orbits with a random (thermal) J-distribution. There is to date 
no satisfactory explanation for the presence of the S-stars so close to the MBH (see Alexander 
|3] for a review). 

The existence of coherent dynamical structures in the GC constrains the relaxation processes 
on these distance scales, since the the relaxation timescales must be longer than the structure 
age t* to avoid randomizing it. Figure (0) shows that the observed systematic trends in the 
spatial distribution, age and state of relaxation of the different stellar components of the GC are 
consistent with, and perhaps even caused by RR. The star disks are young enough to retain their 
structure up to their inner edge at 0.04 pc, where t* ~ -^~Ptp{, vector RR can randomize the disk 
(Hopman & Alexander |27j). It is tempting to explain the S-stars as originally being the inner 
part of the same disks that are currently present in the GC. However, this scenario is somewhat 
problematic. First, we note that vector relaxation can only change the inclinations of the orbits, 
and not their eccentricities, while many of the S-stars have high (e > 0.9) eccentricities; the 
scalar resonant relaxation time is larger than the age of the disks. Second, resonant relaxation 
alone cannot explain why the S-stars are systematically less massive than the disk stars. An 
alternative (Levin [23]) would be that the S-stars were perhaps formed in previous accretion 
disks of which the dynamical signatures have now disappeared. 

If the S-stars were not formed in the disk, but captured by either a tidal binary disruption 
(Gould & Quillen [22]; see also contribution from Perets et al. in this volume) or an exchange 
interaction with a stellar mass black hole (Alexander & Livio [3]), they may be much older than 
the disks, and in particular their age may be comparable to the local scalar RR time (see figure 
In this case, RR will redistribute their orbits within their life-time. This may be an essential 
element of these formation mechanisms: both scenarios lead to rather eccentric orbits (especially 
tidal binary disruption), whereas not all the orbits of the S-stars are very eccentric: star SI has 
eccentricity e = 0.358 ± 0.036, and S31 has e = 0.395 ± 0.032 (Eisenhauer et al. ^3]). Since the 
age of these stars may well exceed the RR time, RR may have redistributed the eccentricities 
to the current DF, which is consistent with a thermal DF. 

Regardless of the origin of the S-stars, their random orbits are consistent with the effect of 
RR. Vector RR can also explain why the evolved red giants beyond 0.04 pc, in particular the 
more massive ones with t* <Cmin(TNR, ? RR ) are relaxed, since T RR <i* out to ~1 pc. 

3. Gravitational wave sources 

MBHs with masses M, < 5 x 1O 6 M have Schwarzschild radii r$ = 2GM,/c 2 , such that a test 
mass orbiting at a few r$ emits gravitational waves (GWs) with frequencies 10~ 4 Hz < u < 1Hz, 
detectable by the planned space based Laser Interferometer Space Antenna? (LISA). Such GW 
sources, for which the mass of the inspiraling object is many orders of magnitude smaller than 
the mass of the MBH are known as extreme mass ratio inspiral sources (EMRIs). GW inspiral 
events are very rare (of the order of 10" 7 - 10~ 8 yr 1 per galactic nucleus; e.g. Hils & Bender |23j : 

3 http://lisa.jpl.nasa.gov/ 



Sigurdsson & Rees (SOj ; Ivanov [22] ; Freitag JZj ; Alexander & Hopman [2] ; Hopman & Alexander 
[26] . [27| . [28]). and it is unlikely that we will observe GWs from our own Galactic center (GC), 
although such a possibility is not entirely excluded (Freitag 18 ; Rubbo, Holley-Bockelmann 
& Finn @Ej)- The galactic MBH plays nevertheless a role of importance in understanding the 
dynamics of EMRIs, since its mass is very close to the mass of the "optimal" LISA EMRI target, 
and as a consequence one may use the GC to model extra- galactic nuclei. 

Hopman & Alexander j^H] used a model based on the GC to analyze the dynamics of EMRIs. 
One of the main results was that inspiraling stars always originate very near the MBH, within 
a distance of ~ 0.01 pc: due to the relatively short relaxation time in galactic nuclei, stars that 
start to spiral in from larger distances are very likely to plunge into the MBH before becoming 
observable as GW emitters. This result was confirmed qualitatively by iV-body simulations by 
Baumgardt et al. 9,; of tidal capture of MS stars by an intermediate mass black hole (Hopman, 
Portegies Zwart & Alexander [21]; Hopman &; Portegies Zwart |25|). 

The fact that only stars within ~ 0.01 pc spiral in successfully, implies that it is the stellar 
content and dynamics of that region which determine the rate of GW inspiral events for the 
different populations in the system. This means, for example, that mass-segregation is likely to 
play an important role (Hopman & Alexander |^Hj; Freitag et al.; |19j : see also contribution by 
Marc Freitag in this volume). Since the resonant relaxation time is very short (Trr <C Tnr) 
near ~ 0.01 pc, it also implies that RR will dictate the rate at which stars are driven towards 
low J orbits, where energy dissipation is efficient and stars spiral in. 

Hopman & Alexander [2Z] used a Fokker-Planck method in energy space with a sink term due 
to RR losses in J-space to calculate the GW inspiral rate. At every time-step, stars redistribute 
in energy-space due to (non-resonant) two body scattering, and stars are accreted by the MBH 
with some specified rate per energy bin. In the relevant regime, this rate is assumed to be 
of order ~ N(E)/Trr, i.e., within one RR time all stars in the bin would be accreted if they 
were not replaced by new stars that flow to higher energies (tighter orbits). In spite of the fact 
that stars are drained very efficiently near the MBH, Hopman & Alexander |22j found that the 
rate at which stars are replenished by two-body scattering is sufficiently high that the stellar 
distribution will not be depleted near the MBH, unless the efficiency of RR is more than an 
order of magnitude larger than exploratory ./V-body simulations (Rauch &; Tremaine [1H]) have 
indicated. Modifications of the stellar DF due to RR are too small to be observable. Some of the 
stars that are captured by the MBH are swallowed directly without giving a GW signal, but the 
stars closed to the MBH (within ~ 0.01 pc) will spiral in rapidly enough to give obtain an orbit 
of period P < 10 4 s for more than a year. Such sources would be observable to LISA to distances 
up to a few Gpc, depending on the mass of the inspiraling star. Since the enhanced rate at which 
stars flow to the loss-cone in angular momentum space is sustained by the larger flow due to 
two-body scattering in energy space, the rate at which EMRIs are produced is increased. The 
analysis by Hopman & Alexander [27] indicates that the rate at which observable EMRIs are 
formed in galactic nuclei is ~ 8 times higher than that for the case in which RR was neglected. 

4. Conclusions 

Resonant relaxation is a relatively unexplored dynamical mechanism. Vector RR, which only 
affects the orientation of the orbit but not the eccentricity, operates in many stellar systems, 
while scalar RR, which does affect the eccentricity, is unique for MBH systems. This mechanism 
becomes important at distances < 0.1 pc from the MBH. It may have played an important role 
in redistributing the orbits of the S-stars, and enhances estimates of the GW inspiral rate by 
nearly an order of magnitude. Our own GC provides a unique case study for resonant relaxation. 
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